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Abstract: The paper proposes a framework for the construction of solutions to a hybrid 
dynamical system that exhibit Zeno behavior. A new approach that enables solution to be 
prolonged after reaching its Zeno time is developed. It allows for a comprehensive stability 
analysis and asymptotic behavior characterization of such solutions. The results are applicable 
to a wide class of hybrid systems and match with practical experience of simulation of real-world 
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1. INTRODUCTION 

Processes that combine continuous and discontinuous be¬ 
havior naturally arise in a variety of real-world applica¬ 
tions such as robotics, biological systems, chemical ki¬ 
netics, logistics and networked control systems. The ba¬ 
sic framework to model and analyse such a behavior 
is impulsive differential equations Samoilenko and Per- 
estyuk (1987); Lakshmikantham et al. (1989); Samoilenko 
and Perestyuk (1995). Besides this theory we would 
like to mention the other more recent developments in 
the related fields like hybrid dynamical systems Van 
Der Schaft and Schumacher (2000), dynamic equations on 
time scales Bohner and Peterson (2012), discontinuous dy¬ 
namical systems Akhmet (2010), switched systems Liber- 
zon (2012) and hybrid automata Henzinger (2000). 

Throughout the paper we will use one of the most recent 
and rapidly developing framework — a hybrid dynamical 
system proposed in Goebel et al. (2012). This framework 
is one of the most general and includes a majority of other 
classes of systems that model processes with continuous 
and discontinuous behavior. Moreover a variety of novel 
results are developed in Goebel et al. (2012) that are not 
available in the other frameworks. Also this framework 
appears well-adapted to the control-related problems. In 
particular the introduction of input-to-state stability (ISS) 
concept for hybrid dynamical systems gave a strong push 
and motivated a fast development of new methods for sta¬ 
bility analysis of hybrid systems with exogenous input Gai 
and Teel (2009). The questions on robustness of ISS for 
hybrid systems were considered in Gai and Teel (2013). 
In recent years a considerable attention is paid to the 
stability analysis of interconnections of hybrid dynamical 
systems. Small-gain approach proved to be an effective 
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tool for stability analysis of solutions to interconnections 
and networks of a large scale Sanfelice (2011); Dashkovskiy 
and Kosmykov (2013); Mironchenko et al. (2014); Sanfelice 
(2014); Liberzon et al. (2014). In spite of these develop¬ 
ments, interconnections of hybrid systems are considered 
only under strong constraints, that are often not compat¬ 
ible with applications Sanfelice (2011); Dashkovskiy et al. 
(2013). 

The simplest example of unsolved problem is related to a 
bouncing ball modelled by hybrid dynamical system. The 
origin of the bouncing ball system is in some sense asymp¬ 
totically stable (for a precise definition see Definition 2.8). 
There is a variety of Lyapunov-like theorems in Goebel 
et al. (2012) to verify this. However if we consider two 
such balls as one system (a so-called vacuous interconnec¬ 
tion) then there are no methods to prove the asymptotic 
stability of the origin for the entire system. Moreover, 
this system is not asymptotically stable in the framework 
of Goebel et al. (2012). It seams that this framework is 
not suitable for modeling this rather simple mechanical 
system, however we claim that the stability problem can 
be resolved by a minor extension of the theory developed in 
Goebel et al. (2012). For more details of the just mentioned 
problem we refer to Section 3 and Figure 1. Here we only 
mention that this problem is caused by the Zeno solutions 
characterized by infinitely many impulsive jumps over a 
finite period of time. Such solutions are not defined after 
this time period. 

Several approaches were proposed to cope with this prob¬ 
lem. Some of these methods enable a solution to be pro¬ 
longed beyond its Zeno time but only for certain classes 
of hybrid systems. In Johansson et al. (1999), a so-called 
regularization technique has been proposed and was illus¬ 
trated for particular examples. It is based on perturbing 
the hybrid system in order to obtain non-Zeno solution, 
and then taking the limit as the perturbation goes to 
zero. A more formal procedure for obtaining generalized 



solutions of Zeno hybrid system via regularization was 
presented in Goebel et al. (2004); Sanfelice et al. (2008). 
For a particular class of Lagrangian hybrid systems, a 
solution switches to a holonomically constrained dynam¬ 
ical system after the Zeno point is reached Ames et al. 
(2006), Or and Ames (2011). In the closely related class of 
switched systems Liberzon (2012), Shorten et al. (2007), 
a solution may converge to a switching surface in a finite 
time, along with increasingly fast switching events near 
this surface. This phenomenon is called chattering. In this 
case, the solutions can be extended by considering the 
set-valued Filippov solution Filippov and Arscott (1988), 
which involves sliding along the switching surface. In Cui- 
jpers et al. (2001), a solution prolongation beyond Zeno 
was proposed by introducing the concept of transition over 
infinite sequence and accumulation-closed transition sys¬ 
tems. Finally, considerable achievements were made from 
a computer science viewpoint. The existence of infinitely 
many discrete events over a finite period of time force 
simulators to ignore some events or looping indefinitely. 
The ways to overcome these problems were proposed in 
Konen et al. (2016) by introducing new algorithms for 
event detection and localization. 

A peculiarity of hybrid dynamical systems is that the 
concept of time is characterized by two parameters: the 
amount of time passed and the number of jumps that 
have occurred. According to Goebel et al. (2012), a certain 
subset of M>o x No is called hybrid time domain. More 
general rules for constructing hybrid time domains were 
proposed in Gollins (2006); Davoren and Epstein (2008). 
In Collins (2006), the concept of generalized hybrid time 
domain has been introduced where a discrete-time axis was 
generalised to a countable ordinal that can have infinitely 
many accumulation points, which correspond to Zeno 
occurrences. This approach enables to prolong solutions 
to a hybrid system beyond Zeno time. However in Collins 
(2006); Davoren and Epstein (2008) authors did not study 
stability properties of prolonged solutions which is in the 
main focus of our paper. 

The aim of this paper is to develop an approach for 
solutions prolongation over the Zeno time and to study 
their stability properties. In view of the well developed 
stability theory in Goebel et al. (2012) we aim to introduce 
some minor extensions in its framework so that we still can 
use results from Goebel et al. (2012). For this reason we 
do not follow such deep modifications of the hybrid time 
domain notion as in Collins (2006); Davoren and Epstein 
(2008). Our slight extension enables stability analysis of 
hybrid systems beyond Zeno points. 

In this paper we propose an approach to extend a solution 
to a hybrid dynamical system beyond its Zeno time with¬ 
out destroying the key concepts of Goebel et al. (2012). 
In our mind, a natural way is to prolong Zeno solution 
from its w-limit point. For this purpose we adapt hybrid 
framework from Goebel et al. (2012) by introducing a 
three dimensional hybrid time domain and redefining the 
concept of solution. 

The rest of the paper is organized as follows. In Section 
2 we recall some basic definitions from the theory of 
hybrid dynamical systems. A motivating example is given 
in Section 3. A new approach for solution construction 


is presented in Section 4. In Section 5 we prove a series 
of propositions that enable stability analysis of solutions 
to a hybrid dynamical system with Zeno behavior. An 
illustrative example is given there. A short discussion on 
open problems in Section 6 completes the paper. 

2. PRELIMINARY NOTION AND DEFINITIONS 

The following notation and definitions are taken from 
Goebel et al. (2012): 

i = f{x), xeC, 
x~^ = g{x), X G D. 

The state a; G M", n G N can change according to the 
differential equation x = f{x) while a; G C, and it can 
change according to the difference equation x'^ = g{x) 
while X € D. The sets C C M" and D C M” are called the 
flow and the jumps sets respectively, functions / : G —)■ M" 
and g : D —>■ M" are the flow and jump maps. The data of 
the hybrid system H is given by {C, f,D,g). 

The parametrization of a solution to the hybrid system 
H is given by two parameters: t G K>o = [0,oo), the 
amount of time passed, and j G No = N U{0}, the number 
of jumps that have occurred. A certain subset of M>o x Nq 
can correspond to evolutions of hybrid systems. Such sets 
are called hybrid time domains. 

Definition 2.1. (Hybrid time domain). Let <ti <t 2 < 
to < .... A subset 

E = C M>o X No 

3 

is a hybrid time domain if it is a union of a finite or infinite 
sequence of intervals [tj,tj+i] x {j}, with the last interval 
(if existent) possibly of the form [tj,T) with T finite or 
T = oo. 

Given a hybrid time domain E we denote: 

sup A = sup{t G IR>o : 3j G No such that (t, j) G E}, 

t 

sup A = sup{j G No : 3t G M>o such that (t, j) G E}. 

3 

Definition 2.2. (Hybrid arc). A function ^ : A —>• M" is 
a hybrid arc if A is a hybrid time domain and if for 
each j G N, the function t —)• (j){t,j) is locally absolutely 
continuous on the interval N = {t : {t,j) G E}. 

Given a hybrid arc </>, the notation dom^ represents its 
domain, which is a hybrid time domain. 

Definition 2.3. (Complete hybrid arc). A hybrid arc f : 
E —>■ M" is called complete if dom0 is unbounded, i.e., 
if supj E + supj E = oo. 

Definition 2.f. (Zeno hybrid arc). A hybrid arc f) : E ^ 
M” is called Zeno if it is complete and sup^ domf) < oo. 

The existence of a Zeno hybrid arc means that an infinite 
number of jumps occurs during a finite time. The time 
T = sup( dom (j) is called a Zeno time. 

Definition 2.5. (Solution to a hybrid system). A hybrid arc 
(j) is a solution to the hybrid system H if 0) G C U D 
and 

(SI) for all j G N such that D := {t : {t,j) G dom(()} has 
nonempty interior 



4>(t,j)&C for all fGintP, 

4>{tjj) = for almost all t G P; 

(S2) for all {t,j) G dom^ such that {t,j + 1) G dom(/), 

G D, (pitj + 1) = g{(t){t,j)). 

The properties of hybrid arcs (like completeness, Zeno, 
etc.) are automatically extended on the corresponding 
solutions. 

Definition 2.6. (Maximal solution). A solution ()) to H is 
maximal if there does not exist another solution ip io % 
such that dom^ is a proper subset of domfi and = 

for all {t,j) G dom^. 

Let S-uiA) denote the set of all maximal solutions </> to a 
hybrid system H with ^(0, 0) € A. 

Definition 2.1. (Strong forward pre-invariance). A set A C 
M" is said to be strongly forward pre-invariant (SFpI) if 
for every p G 5^(A), rge p C A, where rge p = {y gM.^ : 
3(t,j) G dom(() such that y = (j){t,j)}. 

For a precise definition of stability we recall the definitions 
of standard functions and distance to a closed set. A 
function a : M>o -G M>o is called a class-/Coo function (a G 
/Coo) if 01 is zero at zero, continuous, strictly increasing, and 
unbounded. A function p : IR>o -G M>o is positive definite 
{p G W) if p{s) > 0 for all s > 0 and /9(0) = 0. Given a 
vector a: G M" and a closed set A C M", the distance of x 

to A is defined by IxU := inf \x — y\. 

veA 

Definition 2.8. (Uniform global pre-asymptotic stability). 
Let A C M" be closed. The set A is said to be 

• uniformly globally stable (UGS) if there exists a 
function a G /Coo such that any solution p to H 
satisfies \4>{t,j)\A ^ cr(l^^(0)0)U) for all {t,j) G 
domp; 

• uniformly globally pre-attractive (UGpA) if for each 
e > 0 and r > 0 there exists T > 0 such that, for any 
solution (ptoH with |0(O, 0)|^ < r, (/, j) G dom^ and 
t + j>T imply |0(/, j)U < e; 

• uniformly globally pre-asymptotically stable (UG- 
pAS) if it is both uniformly globally stable and uni¬ 
formly globally attractive. 

Definition 2.9. (w-limit set of a hybrid arc). The w-limit 
set of a hybrid arc p : damp -A M", denoted D,{p), is the 
set of all points x G M” for which there exists a sequence 
{(/, j)i}“i of points (ti,ji) G dom^ with lim fi + ji = oo 

i—¥oo 

and lim p{ti,ji) = x. Every such point is an w-limit point 

i^oo 

of p. 


3. MOTIVATING EXAMPLE 

Gonsider two hybrid dynamical systems Hi with states 
Xi G M"‘ and inputs ui Gp C IR™’ 

G Gj, (fi~L P 

where ni,mi G N, i = 1,2. The sets Ci C M”* x Ui 
and Di C M”* x Ui define the flow and the jumps sets 
respectively, functions fi ■ Ci M"* and pi : Di ^ K"’ 
are the flow and jump maps. The data of the hybrid system 
Hi is given by (Ci, fi, D^, gfi. 


Let us interconnect these two systems with ui = hi(xi) 
and U 2 = h 2 {x 2 ), where functions hi : -A U 2 , 

/12 : —>■ Ki. Then the entire interconnection can 

be represented as a single hybrid dynamical system H 
with data (C, f,D,q), where its state is x := (xi,X2) G 
IR”i X IR”U its flow set is 

C :={x: {xi,h2{x2)) G Ci} O {x : {x2, hi{xi)) G G 2 }, 

its flow map is f{x) := {fi{xi,h2{x2)), f2ix2,hi{xi))), its 
jump set is 

D := {x : (a;i, h2{x2)) G Di} U {x : {x2, hi{xi)) G £> 2 } 

and its jump map is 6 f(a;) := (gi(a:i,/i 2 (a; 2 )), 52 ( 2 : 2 ,/ii(a;i))) 
with 

~ \ \9i{xi,h2{x2)), if {xi,h2{x2)) G Di, 

^ otherwise. 


(92{x2,hiixi)), if ix2,hiixi)) G D 2 , 
( 2:2 otherwise. 


In the literature Dashkovskiy and Kosmykov (2013); 
Dashkovskiy et al. (2013) such choice of the flow set C 
and the jump set D is called natural. An important fact 
is that an interconnection of two hybrid systems "Hi and 
H 2 is a hybrid system of the form %. So one may use a 
variety of previously developed methods and techniques 
(for instance from Goebel et al. (2012)) for a qualitative 
characterization of solutions and the problem of a compre¬ 
hensive analysis of interconnections seems to be solved. 
However an essential problem appears in this context. 
It was discussed in Sanfelice (2011) and caused by the 
interconnection of a hybrid system with Zeno solution and 
a hybrid system with continuous complete solution. Such 
interconnection has a Zeno solution that is not a part 
of the set of solutions to every subsystem. Another good 
illustration of this problem is a vacuous interconnection of 
several bouncing balls when the balls start from different 
initial positions. The solution of such model may not allow 
all the balls to reach their own Zeno time as the original 
model of each bouncing ball does (see Figure I). This leads 
to unnatural loss of asymptotic stability of the origin. 

In this paper we propose a way to extend the hybrid 
framework Goebel et al. (2012) in order to cope with 
aforementioned problems. 


4. HYBRID FRAMEWORK EXTENSION 


The main source of the problems stated in the motivation 
section is that a solution to a hybrid system is not defined 
beyond its Zeno time. However some experiments from 
real life like bouncing ball argue that a solution should 
be prolonged over its Zeno time. A bouncing ball after 
reaching the resting state continues to lie while time is 
counting further and further. This motivates us to allow 
solution to continue its evolution after reaching Zeno time. 
In our extended framework, Zeno solution continues its 
evolution from an w-limit point after reaching its Zeno 
time. It enables us to construct solutions that reflect real- 
world observations and to perform their stability analysis. 

To describe the evolution of solution to a hybrid system we 
introduce a new notion of hybrid time domain. It tracks 
not only the elapsed time and the number of impulsive 
jumps, but also the number of Zeno points occurred 




Fig. 1. Evolution of the height coordinates (top) and 
hybrid time domain (bottom) of two vacuously in¬ 
terconnected bouncing balls started with the height 
3 (red) and 1 (blue) respectively. The state of the 
system converges to some point away from the origin 
when hybrid time goes to infinity {t + j —>■ oo). The 
solution is not defined beyond the observed Zeno time. 

during the evolution process. Similar to a classical hybrid 
time domain from Definition 2.1, only certain subsets 
of K>o X No X No can correspond to evolutions of hybrid 
systems. 

Definition fil. (Extended hybrid time domain). Let 
be a set of time moments such that tj^k < tj+i,k Vj, k gNq 
and tj^k < ti,fc+i Vi, j, k G No. A subset 

^ C K>o X No X No 

j,k 

is an extended hybrid time domain if it is a union of a finite 
or infinite set of intervals [tj^k, tj+i,k] x {i} x {k}, with the 
last interval (if existent) possibly of the form [tj^k, T) with 
T finite or T = oo. 

Index k corresponds to the number of encountered Zeno 
behaviors. For a given extended hybrid time domain E we 
denote: 

sup E = supjfc G No : G ffi>o, j G No s.t. {t,j,k) G E}. 

Zeno 

Note that for any extended hybrid time domain we can 
fix an admissible index k and consider ist subset corre¬ 
sponding to this k. Its projection onto K>o x No (dehned 
by dropping k) is the ’’classical” hybrid time domain from 
Definition 2.1. 

An extended solution 0 is a function defined on an extended 
hybrid time domain. Before reaching the first Zeno time 
the extended solution (j) coincides with the ’’classical” 
solution (j) to a hybrid system: (p{t,j,0) = fi{t,j) for all 
(t,j) G dome/). 

In the original framework Goebel et al. (2012) a state 
x{t) = <p{t,j) G M" can evolve along a trajectory of 
differential equation x = f{x) while x(t) G C. At the time 
s G M>o when x(s) G D it can be instantly transferred into 
a new position g{x{s)) and the value of the corresponding 
jump index in hybrid time domain increases by 1 so that 


ip{s,j -I- 1) = g(x(s)). In our settings we add one more rule 
to construct solution to a hybrid dynamical system: 

• if for some fixed fc G No hybrid arc (j){t,j,k) is Zeno 
with non-empty w-limit set then a solution to a 
hybrid system % is prolonged with initial condition 
<j>{T,0,k -I- 1) G fl(0), where r is the Zeno time for 
hybrid arc (/>(•,•, fc). 

Our extended solution is a concatenation of classical 
hybrid arcs (j)^{t,j) with initial conditions ((>°(0,0) = 

G ^^(r2,0) G and so on, where is 

the Zeno time for the hybrid arc . 

A new rule of extended solution’s construction leads to the 
following properties of the corresponding extended hybrid 
time domain E : if the point (to,fe-i-ij Oj ^+1) G E then there 
exist infinitely many points of the form (t.^fc, •, fc) G E such 
that lim tj^k = to,fc+i- 

j-foo 


In general, an w-limit set D.{(j)) may consist of several 
or infinitely many points. According to our new rule a 
single initial point can generate multiple solutions. Such 
situation appears, for example, in modelling of water tanks 
system (see Alur and Henzinger (1997) for details). This 
system has Zeno arcs with two w-limit points. Therefore 
two different extended solutions will be generated from 
a single initial point. This is quite natural since the 
considered physical process can evolve according to both 
of solutions in a real experiment. 

Remark J^.l. For a given Zeno hybrid arc ip the w-limit 
point X G fl(¥’) is a limit of a sequence of points {xi\ of 
the state space such that Xi G D, i = 1,2,.... A behavior of 
the corresponding extended solutions now heavily depends 
on the properties of the jump set D. If D is a closed 
set (it means that it contains all its limit points) then 
the extended solution continue its evolution from limit 
point X G D and therefore should jump. This can lead to 
eventually discrete solution. If the jump map D is an open 
set then the extended solution can continue its evolution 
from the limit point x G C and therefore can be prolonged 
beyond Zeno time of the ordinary time axis. From this 
viewpoint it is quite natural to model real processes with 
an open jump set D as it is done in the following example. 

Example fil. Consider a vacuous interconnection of two 
bouncing balls. Let xi,X3 G K>o stand for the heights 
of the balls and X2,X4 G K stand for the corresponding 
velocities. Then system has the form 

Xi X 2 , X 3 X 4 , X G C 

X2 = -7(xi,X2), X4 = -7(x3,X4), 


Xi+ = Xi, 


= X 3 , 


X2~^ = 


' —Xx 2 ,x G Di, 


\X 2 , 


X ^ Di 

X4 

\X 4 , 

X ^ D2, 

Cl = 

{x 

G K"* : 

Xi > 

II 

u 

0 

0 

0 ,X 2 > 0 }, 

Di = 

{x 

G : 

Xi = 

0 ,X 2 < 0 } 


C 2 = 

{x 

G K"* : 

X3 > 

0 or X3 = 

0, X4 > 0}, 

D2 = 

{x 

G : 

X3 = 

0, X4 < 0} 

5 

C = 

Cl 

nC2, 

D = 

= Di U £>2 









where A G (0,1) is the restitution coefficient, 7 : —>■ K 

is given by 



Fig. 2. Evolution of the height coordinates (top) and ex¬ 
tended hybrid time domain (bottom) of two vacuously 
interconnected bouncing balls started with the height 
3 (red) and 1 (blue) respectively. 

A numerical simulation is presented on Figure 2. The arc 
that corresponds to Zeno index k = 0 fully coincides with 
the one from the original framework Goebel et al. (2012). 
Its cu-limit set consists of a single point (in this case the 
uniqueness of solution is preserved). At the Zeno time of 
the blue ball the solution is now prolonged from this point 
and Zeno index is increased by 1. The extended solution 
exhibits a further Zeno behavior and its w-limit set is just 
the origin. At the Zeno time of the red ball, the solution 
is prolonged from its w-limit set (0,0,0, 0) S which is 
a single point again. The last arc of solution is trivial and 
purely continuous with supj dom^ = 00. The concatenated 
solution corresponds to our experience. 

5. STABILITY ANALYSIS 

In this section we introduce a new stability notion in order 
to describe asymptotic behavior of extended solutions to 
hybrid systems. Two auxiliary lemmas will be needed to 
justify stability characterization. 

Lemma 5.1. UGS of a set A implies its strong forward 
pre-invariance. 

Proof. Suppose it is not true. Let there exist a solution 
(p to Ti. with 0(0,0) S A and a point x* G rge0(t,j) 
such that X* ^ A for some (t,j) ^ dom0. It means that 
there exists <5 > 0 such that |0(t,J)U = 5. Then from 
Definition 2.8 it follows directly that there is a function 
a G /Coo such that 

0 < d = |0(t, j)U < a(l0(O, 0)U) = a(0) = 0. 

The contradiction proves that every solution 0 starting in 
A remains in this set: (j){t,j) G A for all (t,j) G dom0. 
Lemma 5.2. Let A be UGpAS and every arc 0 with initial 
condition in {C'UZ/}\Al have a non-empty w-limit set, 
then n(0) C A. 


Proof. Suppose it is not true. Consider a solution 0* to 
Ti. with 10* (0,0)1^ < r that has an w-limit point ^ outside 
the set A. It means that there exists a sequence {(t, j)i}^i 
of points G dom0 with lim 0 + = 00 and 

i—^oc 

lim 0*(0,ji) = C ^ A. (2) 

Then there exists d > 0 such that = S. From the 
UGpAS of the set A it follows that for e = | there exists 
T > 0 such that for every (t,j) G dom0 with t + j > T 
follows 10(0 j)U < 5- 

However the existence of the limit (2) guarantees that 
for any d > 0 there exist a (/-neighbourhood Ud{0 and 
G dom0 such that 0(/*,j*) e Ud{^). Choosing 
d small enough to satisfy the conditions t* + j* > T 
and Ud{^) n Us (A) = 0 leads to |0(/*, j*)U > 2’ which 
contradicts the UGpAS of the set A. This proves that the 
w-limit point ^ G A. 

Definition 5.1. ("Hn A). If A C C U D is UGpAS for %, 
then A can be considered as the state space for a new 
hybrid system with the new flow set C n A and the new 
jump set Df^A. We will denote this new system by F/ H A. 

Indeed, from Lemma 5.1, UGpAS implies SFpI of the set 
A so every solution with initial condition in A will remain 
there for all (/, j) G dom0. In the case when 0 is a Zeno 
solution it will be prolonged from a point of its w-limit 
set H(0). From Lemma 5.2 follows that H(0) C A, so 
the solution will again remain in the set A. It means that 
extended solutions of the system 7/ with initial conditions 
in A will coincide with extended solutions of the system 
7/ n A with the corresponding initial conditions. Since 
A C C U 77, no new solution will be generated. 

For a comprehensive description of asymptotic behavior 
of extended solutions we introduce a new definition of 
stability over Zeno. 

Definition 5.2. (UGpASoZ). Let A S K" be closed. The 
set A is said to be 

(i) uniformly globally stable over Zeno (UGSoZ) if there 
exists a function a G /Coo such that any solution 0 
to 7/ satisfies |0(/,7,/c)U ^ Q^(l0(Oi 0i 0)U) for all 

(t,J,k) G dom0; 

(ii) uniformly globally pre-attractive over Zeno (UG- 
pAoZ) if for each e > 0 and r > 0 there exist 
T > 0 and K > 0 such that, for any solution 0 to 
7/ with 10(0,0,0)1^ < r, from {t,j,k) G dom0 with 
either t + j >T, k = K 01 k>K or t + j>T, 
k = supzeno dom0, sup^eno dom 0 < 7C it follows that 

I0(/, j. fc)U < e; 

(iii) globally pre-attractive over Zeno (GpAoZ) if for each 
e > 0, r > 0, and for any solution 0 to 7/ with 
10(0,0,0)1^ < r, there exist T > 0 and K > 0 such 
that from (/, j, k) G dom0 with either t+j > T, k = K 
or fc > 7F it follows that |0(/, j, fc)U < e; 

(iv) uniformly globally pre-asymptotically stable over 
Zeno (UGpASoZ) if it is both UGSoZ and UGpAoZ. 

The conditions for the pre-attractivity actually mean that 
all solutions will reach the e-neighbourhood of the set A 
no later than at the time T after K-t\i Zeno occurrence. 
The uniformity means that T and K are the same for all 









solutions. If a solution does not undergo such number (K) 
of Zeno occurrences then it should reach the corresponding 
e-neighbourhood no later than time T after its last Zeno. 

Theorem 1. Let there exist a finite sequence An C An-i C 
... C Ai C Aq = C yj D such that Ai is UGpAS for the 
system TiCiAi-i, i = and for all initial values 

^(0,0) € Ai-i\Ai, i = l,...,n — 1 the corresponding 
solutions (j) are Zeno with non-empty w-limit sets. Then 
An is UGpASoZ. 

Proof. If n = 1 then UGpAS of the set Ai implies its 
UGpASoZ with AT = 0. Let us consider the case n = 2. 
First we prove stability of the set A 2 - From UGS of 
the sets Ai and A 2 it follows that there exist functions 
01,02 € Aoo such that 

< oi(|^(0,0)UJ Vv3(0,0) G Ao, 
\‘pit,j)\A2 < a2(|«i'(0,0)UJ Vv3(0,0) G Ai 

and for all {t,j) G dom^. Then the extended solution (f 
satisfies 

\Tit,j,k)\A2 < ai2(|^(0,0,0)UJ V(^(0,0) G Aq 

for all G dome/) with 012 ( 5 ) = max{oi(s), 02 ( 5 )}. 

Hence A 2 is tlGSoZ. 

Now we will prove the pre-attractivity. For this purpose 
we will show that each extended solution (j) to the hybrid 
system TL satisfies the uniform pre-attractivity conditions 
{ii) of Definition 5.2 with respect to the set A 2 . Note that 
since A 2 is UGpAS for the system "HnAi every solution 
initiated from Ai satisfies the pre-attractivity conditions 
of the Definition 2.8: Ve,r > 0 there exists Ti > 0 such 
that any solution to "HnAi such that |^(0, 0)|^2 < r, 
^(0,0) G Ai satisfies \(l){t,j)\A 2 — ^ ^ dome/) 

with t + j > Ti. It means that the corresponding extended 
solution satisfies the uniform pre-attractivity conditions 
(ii) of Definition 5.2 with T = Ti and K = 0. 

It remains to show that the conditions (ii) of Definition 5.2 
are also satisfied for solutions starting outside the set 
Ai- Let |(^(0,0,0)1^2 < r, ^(0,0,0) G Ao\Ai and let 
the arc (f) = € dom(p be Zeno. From 

the conditions of Theorem 1 its w-limit set is non-empty, 
hence this solution is being prolonged from the set 
From Lemma 5.2 it follows that C Ai and from 

Definition 5.1 it follows that the set Ai can be considered 
as a new state space for the system 'HCiAi. Since A 2 is 
UGpAS for the system HCiAi and D((/)) C Ai it follows 
that an extended solution (j) issued from Ao \ Ai satisfies 
the pre-attractivity conditions (ii) of Definition 5.2 with 
T = Ti an a: = 1. 

Since there are no other types of solutions to the system TL 
starting outside Ai, the set A 2 is UGpAoZ with T = Ti, 
AT = 1. UGpASoZ follows from UGSoZ and UGpAoZ. 
Iterating the previous reasoning one can prove UGpASoZ 
for any finite n. This concludes the proof. 

The proven result is applicable only to a system with the 
arcs, issued outside the set Ai, i = 1,... ,n — 1, that are 
Zeno. However it is easy to prove the stability and pre- 
attractivity for the case of non-Zeno hybrid arcs that reach 
the corresponding set Ai in a finite time (e.g. when t + j 
is bounded). In this case we lose the uniformity of the 
pre-attractivity which means that the constants T and 


K that describe the time needed for a solution to reach 
the £-neighbourhood of the set An depend on a particular 
solution. 

Theorem 2 . Let there exist a finite sequence A„ C A„_i C 
... C Ai C Aq = C \J D such that Ai is UGpAS for the 
system TLC\Ai-i, i = l,...,n and for all initial values 
<()( 0 , 0 ) G Ai-i\Ai, i = 1,.. . ,n — 1 the corresponding 
solutions (j) are either Zeno with non-empty w-limit sets or 
rge (/) n Ai 7 ^ 0. Then A„ is UGSoZ and GpAoZ. 

Proof. The proof repeats the reasonings of the Theo¬ 
rem 1. UGSoZ proof is similar. However, we need to 
consider the second type of arcs (non-Zeno) in order to 
prove the pre-attractivity. Let (j) be the solution issued 
from a point outside the set Ai such that rge (() n Ai 7 ^ 0. 
It means that there exists {t*,j*) G dome/) such that 
G Ai. The further evolution of solution 
(j) coincides with the evolution of the solution (p* to the 
system "HnAi with initial condition (/)*(0, 0) = ^*. Since 
the set A 2 is UGpAS for the system "HnAi it follows 
that the solution issued from ^ G Ao \ Ai satisfies the 
conditions (m) from Definition 5.2 with T = t* + j* +Ti, 
K = 0 . Note that t* and j* can be different for every 
solution (j). Denote T2(ip) = t* + j*. 

Since there are no other types of solutions to the system 
TL issued outside Ai, the set A 2 is GpAoZ with T = Ti + 
T 2 ((()), AT = 1. Iterating the previous reasoning one can 
prove GpAoZ for any finite n. This concludes the proof. 

Note that the situation when every set Ai is UGpAS but 
there exists an arc that wounds by spiral and tends to 
the set Ai but does not intersect it, does not satisfy the 
conditions of Theorem 2. 

Next we present an example that demonstrates the usage 
of Theorem 1. To check the UGpAS of a set A we will use 
the known theorem from Goebel et al. (2012): 

Proposition 5.1. (Goebel et al. (2012)). Let A C M” be 
closed. If V{x) is a Lyapunov function candidate for TL 
and there exist ai,a2 G /Coo, and a continuous p G W 
such that 

ai(|a:U) < V{x) < Q:2(|a;U) Vcc G CUDUgiD) 
{VV{x),f{x)) <-p(|a:U) Vx G C 
y{9{x)) -v{x) < -p(\x\a) Vx G D 
then A is UGpAS. 

Example 5.1. Gonsider the following system with state 

X = {xi,X2, X3) G 

Xi = X2, Xi^ = Xi, 

±2 = —'^{xi,X2), X e C, X2~^ = —Xx2, X G D (3) 
is = -Xs, 3^3^ = -a^3, 

and with flow and jumps sets given by 

G = {a; G : xi > 0 or xi = 0, X 2 > 0}, 

D = {a; G : xi = 0, a ;2 < 0}, 

where A G (0,1) and 7 : —)• M is given by (1). 

Let us prove that the origin is UGpASoZ. As one may 
notice, flow and jump sets of system (3) do not depend on 
X3. This system can be interpreted as an interconnection of 
a bouncing ball with state (a;i, a; 2 ) and some other process 
with state X3. Each time when the ball bounces at the floor 
the state variable X3 changes its sign. 



Let function V be defined by 

V{x) = (1 + 0 arctana; 2 ) 

with 

1 _ y 2 

7 = 9,81. 

7r(l + A"') 

The set Ali = {(a; G : xi = 2:2 = 0)} is UGpAS since 
V satisfies Proposition 5.1 with respect to the origin for a 
single bouncing ball Goebel et al. (2012) and the distance 
from a point (ii, 0 : 2 , 2 : 3 ) G to the set Ali coincides 
with the Euclidean norm H-H of the corresponding vector 
(a;i,a; 2 ) G 

|(a;i,X2,X3)Ui = ||(xi,X2)|| = xj + x^. 

Then we arrive to a system of the form (3) with the new 
state space Ai, the new flow set C n Ai = {a; G Ali : 
a;i > 0 U a;i = 0, a :2 > 0} = All and the new jump set 
H n All := {a; G All : a;i = 0 ,a ;2 < 0} = 0. This system 
is purely continuous and the origin (0,0,0) is UGpAS. 
Since all the arcs issued outside the set Ali are Zeno, from 
Theorem 1 it follows that the origin is UGpASoZ. □ 

Theorem 1 proposes a sequential narrowing of the state 
space of a hybrid system. For the last example this process 
can be described with the following sequence of sets: 


^ 7 ^ \ 

{0}x{0}xIR origin yA 


Theorem 1 has been used to prove UGpASoZ of the origin 
without constructing solutions to the system (3) and their 
prolongation explicitly. One may check that for a particu¬ 
lar initial data the corresponding solutions to the system 
(3) have two w-limit points. If ^(0, 0, 0) = (a, b, c) G 
with \/a^ + b^ 7 ^ 0, c ^ 0 then system (3) has Zeno hybrid 
arc with Zeno time t = j (^b + \/b^ + 2ga^ and the 

w-limit set consisting of two points (0, 0, ±c • e~'^). Hence 
the initial point (a, b, c) generates two solutions. Despite 
such complex situation we were able to use Theorem 1 to 
verify UGpASoZ without knowing the exact number of w- 
limit points of Zeno arcs. Moreover, Theorem 1 along with 
Lemma 5.2 can be used for w-limit points localization. 
If one finds a function V satisfying Proposition 5.1 for 
some set Al, then following Lemma 5.2, w-limit points of 
solutions are contained in the set Al. 

6 . DISCUSSION AND OPEN QUESTIONS 

The results presented here are beneficial for construction 
and stability analysis of solutions to hybrid dynamical sys¬ 
tems that exhibit Zeno behavior. The main contributions 
of the paper are the following. First, we have introduced 
the extended hybrid time domain and new prolonged so¬ 
lution concept that heavily relies on the axiomatics and 
notation of Goebel et al. (2012). These extended solutions 
helped us to avoid such undesired effects as freezing of so¬ 
lutions. Second, we propose a generalisation of the attrac- 
tivity concept and prove theorems that provide Lyapunov- 
like sufficient conditions for stability without knowing the 


explicit solution and w-limit points. However, in order to 
apply these theorems one should be able to verify whether 
all hybrid arcs are either Zeno or intersect an appropriate 
set Ali. 

We believe that the proposed way of solution’s prolonga¬ 
tion from its w-limit points can also be achieved without 
the introduction of a new 3-dimensional hybrid time do¬ 
main. However it would cause a significant redefining of the 
basic concepts of hybrid dynamical systems framework. 
An important advantage of the proposed approach is the 
ability to utilize a wide range of previously developed 
results on UGpAS, e.g. from Goebel et al. (2012), for 
stability analysis of extended solutions. 

Several problems have no answers yet and are very exciting 
to be solved. The first one is an extension of the results 
to infinite dimensional setting. This can be described by a 
vacuous interconnection of infinitely many bouncing balls. 
One can easily check that this system has a qualitatively 
different behavior depending on the initial conditions. Let 
the balls are enumerated by index n from I to 00 and each 
ball starts its way with zero velocity and vertical position 
equals to n. Then the Zeno index k of the hybrid time 
domain for this case tends to infinity while the ordinary 
time t —>■ 00 . However if every ball starts with the position 
A then the Zeno index k reaches infinity by a finite 
ordinary time t. If we interconnect each of these systems 
with a purely continuous process that tends to zero (like 
^ = —x) then a solution of the entire interconnection will 
tend to the origin in the first case and will ’’freeze” away 
from the origin in the second one. This situation gives an 
intuition that such kind of systems can be treated using 
some analogues of a local stability concept and needs a 
deeper investigation for a comprehensive analysis of its 
behavior. 

Another challenging issue is an interconnection of a com¬ 
pletely continuous and a completely discrete system. The 
resulting flow and jump sets obtained in ’’natural” manner 
lead to a system with only discrete time domain. The 
examples of such processes are for instance sample-and- 
hold control where a discrete-time algorithm measures the 
state of a continuous time system and updates it. In this 
case an entire interconnected system will have a solution 
with only discrete time and we just lose the continuous 
process. 
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